Solution to the inverse problem for a noisy spherical gravitational wave antenna 



00 
On 
On 



Stephen M. Merkowitz 
INFN Laboratori Nazionali di Frascati, Via Enrico Fermi 40, F00044 Frascati (Roma) ITALY 

(June 10, 1998) 

A spherical gravitational wave antenna is distinct from other types of gravitational wave an- 
tennas in that only a single detector is necessary to determine the direction and polarization of a 
gravitational wave. Zhou and Michelson showed that the inverse problem can be solved using the 
maximum likelihood method if the detector outputs are independent and have normally distributed 
noise with the same variance. This paper presents an analytic solution using only linear algebra that 
is found to produce identical results as the maximum likelihood method but with less computational 
burden. Applications of this solution to gravitational waves in alternative symmetric metric theories 
of gravity and impulsive excitations also are discussed. 



PACS numbers: 04.80.Nn, 95.55.Ym, 04.30.Nk 



I. INTRODUCTION 



Several resonant-mass gravitational wave antennas are now in continuous operation with strain sensitivities of the 
order 10 _21 Hz -1 / 2 |l]]. With further improvements to these detectors and the addition of several large laser interfer- 
ometers now under construction Q , the prospects for gravitational wave astronomy are quite good. The underlying 
non-gravitational physics associated with these detectors is reasonably understood and further improvements can be 
£> ' based on solid technological guidelines. 

Many believe the next generation of resonant-mass antennas will be of spherical shape Q . Confirmed detection 
of gravitational waves will require a coincidence between several detectors, thus the unique features of a sphere may 
play an essential role in a network of gravitational wave antennas. Two important features of a sphere are its equal 
sensitivity to gravitational waves from all directions and polarizations and its ability to determine the directional 
' information and tensorial character of a gravitational wave |Q . 
Q\ . To take full advantage of these capabilities, one needs to be able to interpret the data such a detector will produce. 
Recently, much work has been done to understand the output of a spherical antenna equipped with resonant trans- 
?"r ducers All of these proposals operate on the principle that the response of the transducers can be transformed 

into a quantity that has a one-to-one correspondence with the tensorial components of a gravitational wave. With the 
(^q| measurement of these components, it is possible to solve the inverse problem to obtain the direction and polarization 
• • ■ amplitudes of a gravitational wave. 
. £h ' The soluti 011 to the inverse problem for a noiseless spherical antenna first was outlined in the mid 1970's by Wagoner 
and Paik ||. More recently, the solution for a network of five noiseless bar antennas or interferometers was solved 
by Dhurandhar and Tinto [Q. This method assumed that the detectors were co- located but oriented in different 
directions. This solution is quite elegant because the exact solution can be found using straightforward algebra. Since 
a sphere can be thought of as five bar detectors occupying the same space, this solution can be adapted for a spherical 
antenna ||j9| . In addition, Lobo outlined a procedure that can be used if the correct theory of gravity is not general 
relativity but unknown [|l 0f| . What all these proposals have in common is that they use basic symmetry properties of 
the matrices describing the detectors and their response to a gravitational wave. This makes them intuitive and easy 
to visualize |ll| . 

The solution to the inverse problem in the presence of noise is more complicated. Giirsel and Tinto solved the prob- 
lem for three noisy interferometers using a maximum likelihood method |p_ 2[| . This solution required the measurement 
of the time delay of the signal between widely separated detectors to triangulate the direction of the source. Zhou 
and Michelson showed that the inverse problem for a spherical antenna can be solved in the presence of noise, also 
using a maximum likelihood method j^] . 

What is disappointing about the maximum likelihood method is that the original simplicity of Dhurandhar and 
Tinto's noiseless solution is lost. In addition, an exact solution for the spherical detector was not found, making it 
necessary to solve the problem numerically. This solution can be computationally expensive, especially if the signal- 
to-noise ratio (SNR) is low. For the three interferometer case, the numerical solution often can lead to an incorrect 
estimate at low SNR fll2| . This appears not to be a problem for a spherical antenna; a global maximum usually exists 
and is aligned with the correct direction. This leads us to believe that the problem for a spherical antenna can be 
solved analytically. 
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On experiments with the room-temperature prototype spherical antenna (TIGA) at Louisiana State University, we 
used a procedure to solve the inverse problem for impulsive excitations applied to the sphere surface Jl^,|lJ|. This 
solution was similar to Dhurandhar and Tinto's original solution for gravitational waves, but we used a perturbation 
argument (presented below) to take into account the finite SNR of the experiment. The case of an impulsive excitation 
is more simple than a gravity wave because fewer parameters are involved, however, we show below that this method 
also can be used to find an analytic solution for gravitational waves. 

We begin by reviewing the response of a spherical antenna to gravitational waves in general relativity and show how 
Dhurandhar and Tinto's original method can be applied to solve for the wave direction and polarization amplitudes. 
We then generalize the arguments to any symmetric metric theory of gravity as well as to impulsive excitations. 
In Sec. [I] we show how this technique can be extended to a noisy antenna with independent and equally sensitive 
detector outputs. This solution is found to be equivalent to the maximum likelihood method under the same noise 
requirements. The general approach taken allow this solution to be easily adapted to other types of excitations with 
similar symmetry properties. We conclude the paper with a discussion of the limitations of the solution and possible 
extensions of this method. 



II. DETECTOR RESPONSE OF AN ELASTIC SPHERE 

Dhurandhar and Tinto solved the inverse problem for 5 bar antennas as well as 5 interferometers 0. Others have 
used their method to solve the problem for a spherical antenna |8|^] . Their technique involves constructing a matrix, 
say A, that describes the response of the detector to a gravitational wave. They found that in general relativity the 
eigenvector of A with zero eigenvalue points in the propagation direction of the wave. In the following we will also 
use this concept, but will derive the equations in the context of linear algebra as this will lead us directly into the 
solution for the noisy antenna. For a more complete discussion of the response of an elastic sphere to gravitational 
waves the reader is referred to Refs. [Il5 16 



A. Detector response in general relativity 

A gravitational wave is a traveling time-dependent deviation of the metric perturbation, denoted by H^ v (t). We 
follow a common textbook development for the metric deviation of a gravitational wave, which finds that only the 
spatial components fljj(t) are non-zero, and further can be taken to be transverse and traceless [ fl7|| . This tensor is 
simplified if we initially write it in the "wave-frame," denoted by primed coordinates and indices. This is a coordinate 
frame with origin at the center of mass of the detector and the z' axis aligned to the propagation direction of the 
wave. We restrict ourselves to detectors much smaller than the gravitational wavelength so only the time dependence 
of Hij{t) will have significant physical effects. A general form for the spatial components of the metric deviation in 
the wave-frame can be written as 

H'(t) = 

where h + (t) and h x (t) are the wave amplitudes for the two allowed states of linear polarization and are called the 
plus and cross amplitudes. 

The detector is more easily described in the "lab-frame," denoted by unprimed coordinates and indices, with origin 
at the center of mass of the detector and the z axis aligned with the local vertical. In this frame, the primary physical 
effect of a passing gravitational wave is to produce a time dependent "tidal" force density f GW (x,t) on the material 
with mass density p at coordinate location Xj. This force is related to the metric perturbation by 

f™(x,t) = ±p'£H ij (t)x j . (2) 

3 

It is natural to look for an alternate expression that separates the coordinate dependence into radial and angular 
parts. Because the tensor ffy(i) is traceless, the angular expansion can be done completely with the five second order 
real valued spherical harmonics Yi m <j)), where the index m = 1 . . . 5. We call the resulting time dependent expansion 
coefficients, denoted by h m {t), the "spherical amplitudes" [fig] . They are a complete and orthogonal representation of 
the cartesian metric deviation tensor Hij (t) . They depend only on the two wave- frame amplitudes and the direction 
of propagation. 



h+(t) h x (t) 

h x (t) -h+(t) , (1) 
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To transform the metric perturbation to the lab-frame we perform the appropriate rotations using the y-convention 
of the Euler angles shown in Fig. 0. We denote the rotation about the wave z' axis by a, the rotation about the new 
y axis (7/) by (3, and the rotation about the final lab z axis by 7. The rotation matrix for the y-convention is 



R 



cos 7 cos (3 cos a — sin 7 sin a cos 7 cos (3 sin a + sin 7 cos a — cos 7 sin (3 
— sin 7 cos (3 cos a — cos 7 sin a — sin 7 cos /? sin a + cos 7 cos a sin 7 sin /3 
sin (3 cos a sin /3 sin a cos /? 



(3) 



At this point we arbitrarily set the rotation a about the wave z' axis equal to zero; inclusion of this rotation will only 
"mix" the two polarizations of the wave. The spherical amplitudes can now be written in terms of the polarization 
amplitudes and the source direction 



hi(t) 
h 2 (t) 
h 3 (t) 
hi(t) 
h 5 {t) 



h + (t ) i ( 1 + cos 2 0) cos 27 + h x (i ) cos (3 sin 27 , 
-h+(t)- (l + cos 2 0) sin 27 + h x (t) cos j3 cos 27, 
— h+(t) - sin 2/3 sin 7 + h x (i) sin /3 cos 7, 
sin 2(3 cos 7 + /i x (t) sin (3 sin 7, 



~V3shr ). 



(4a) 
(4b) 
(4c) 
(4d) 
(4e) 



The mechanics of a spherical antenna can be described by ordinary elastic theory. One finds that the eigenfunctions 
of an uncoupled sphere can be written in terms of the spherical harmonics 



^nhn{r,9,(j)) = (a nt (r)r + f3 nt {r)aV)Y em {9,4>). 



(5) 



The radial eigenfunctions a n i(r) and (3 n e{r) determine the motion in the radial and tangential directions respectively 
and depend on the radius a and the material of the sphere j|,[l6) . 

In general relativity, only the 5 quadrupole modes of vibration will strongly couple to the force density of a 
gravitational wave. For an ideal sphere they are all degenerate, having the same eigenfrequency, and are distinguished 
only by their angular dependence. The effective force i 7 i2m(i) that a gravitational wave will exert on a fundamental 
quadrupole mode m of the sphere is given by the overlap integral between the eigenfunctions of the sphere and the 
gravitational tidal force 



F\2m{t) 



*i2m(a:) • f GW (x,t)d 3 x = h m {t)M X a. 



(6) 



Each spherical component of the gravitational field determines uniquely the effective force on the corresponding mode 
of the sphere and they are all identical in magnitude. We can interpret the effective force Fi2 m (t) m each mode as 
the product of: the physical mass of the sphere M, an effective length \a (a fraction of the sphere radius), and the 
gravitational acceleration ^h m (t). The value of the coefficient x depends on the sphere material, but is typically ~ 0.6 

By monitoring the quadrupole modes of the sphere, one has a direct measurement of the effective force on the 
sphere and thus the spherical amplitudes of the gravitational wave. The standard technique for doing so on resonant 
detectors is to position resonant transducers on the surface of the sphere that strongly couple to the quadrupole 
modes. A number of proposals have been made for the type and positions of the transducers |^,hJ|. What all of 
these proposals have in common is that the outputs of the transducers are combined into "mode channels" g m {t) that 
are constructed to have a one-to-one correspondence with the quadrupole modes of the sphere and thus the spherical 
amplitudes of the gravitational wave |l(],[ll| , 



9m (t) oc F 12m (t) oc h m (t). 



(7) 



The mode channels can be collected to form a "detector response" matrix A(t) that in the absence of noise is equal 
to the cartesian strain tensor H(t) in the lab frame 



A(t) 



9i(t) - fas® 
92(t) 

94(t) 



sa(t) 

9s(t) 



94.® 

.93 (i) 

V3-- 



(8) 
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For the remainder of this discussion we drop the notation of time dependence (t) for brevity. 

The strain tensor in the lab frame H is a symmetric traceless matrix. Consequently, it can be orthogonally diago- 
nalized and has an orthonormal set of three eigenvectors. One can construct from the eigenvectors a transformation 
matrix R that diagonalizes H. The matrix R is also orthogonal, thus it can be considered a rotation matrix (it may 
also include a reflection) . The physical interpretation of this transformation is to rotate the lab frame such that the 
z axis points in the direction of the source. The matrix R (the eigenvectors) will tell us the angles of rotation and 
thus the direction of the wave. 

In the wave frame, H' is not normally diagonal but it can be diagonalized by rotating Eq. (Q) about the propagation 
axes using the Euler angle a. a may be a constant or a function of time depending upon the situation. This rotation 
changes the polarization components of the tensor but not the wave direction relative to the lab frame. 

To calculate the rotation matrix R we need to solve the general eigenvalue equation for the strain tensor 

Hx = X x x. (9) 

Since A and H are equal in the absence of noise we are free to substitute A in Eq. (^|) for H . By inspection of Eq. (|l|) 
we see that in general relativity the eigenvector of H with = points in the propagation direction of the wave. 
The direction can be calculated from this eigenvector by recognizing that it corresponds to the last column vector of 
R in Eq. ([|) . Dividing the elements of this column we find 

tan7=-~ (10) 



tan/ 3=y_ — . (11) 
z sin 7 

The unusual minus sign in Eq. (|l^) comes from the use of the y-convention of the Euler angles. Expanding Eq. (|^) 
for \ x = and substituting in a particular choice of matrix elements from Eq. (|§J) we find 

3ff 4 ff3 ~ 2y / 3g 2 g5 noA 



2V3.9135 + 2g? + 3. 9: 
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tan/3 = ± r ^9*9* - 29*92 _1_ (13) 
V 30154 + .9534 + V 3#352 sm 7 

This solution is valid only for a noiseless antenna; it will fail otherwise because we can no longer replace H with A 



and their eigenvectors and eigenvalues will no longer be equal. The ± in Eq. (13) illustrates the unavoidable fact that 
a single sphere cannot distinguish between antipodal sources. This ambiguity is a characteristic of all gravity wave 
detectors, but can be removed by measuring the time delay of the signal between two widely separated antennas. 

Once the direction is calculated we can determine the two polarization amplitudes by taking a linear combination 
of Eqs. (Q). These equations are actually overdetermined so several solutions exist (we have 5 equations but only 4 
unknowns). In the absence of noise any particular solution to them is valid, but in anticipation of the noisy case we 
will take a systematic approach to the solution. 

We need only the angles (3 and 7 to rotate H to H' , so at this point we again set a = 0. The amplitudes are found 
by equating them to the corresponding matrix elements of H' in Eq. ([!]). Again, we may substitute A for H and 
A 1 for H 1 so we have h + = A' n = —A' 22 and h x = A' 12 = A' 21 . A and A' are symmetric so A' 12 and A' 21 will always 
be identical even when noise is introduced. However, no such restriction is placed on A'^ and A' 22 . We will use the 
average (A' n — A' 22 )/2 to calculate h + for reasons that will become clear later. 

Multiplying A' = RF AR for a = and selecting the proper elements we find 

h + = .91 -(1 + cos 2 0) cos 27 — #277(1 + cos 2 (3) sin 27 — g 3 - sin 2(3 sin 7 + 34- sin 2/3 cos 7 + 35 -r- sin 2 j3. (14) 

__ Zi Zi Li 



hx =3i cos (3 sin 27 + g 2 cos (3 cos 27 + g 3 sin (3 cos 7 + 34 sin (3 sin 7. (15) 

These equations can also be derived by taking a linear combination of Eqs. (Q). This is not the only valid solution 
in the noiseless case, but it is particularly symmetric: the coefficients of each component g m is the same as the 
corresponding coefficients of h + or h x in Eqs. (|J) for h m . The fact that h x does not contain a 35 contribution is an 
artifact of using the y-convention of the Euler angles; in other conventions this term may be non-zero. 
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B. Detector response in alternative theories of gravity 



Experiments in the solar system and pulsar-timing tests have ruled out many competing theories of gravity, however, 
general relativity is not the only theory of gravity that passes these weak field tests . One measurement that can 
potentially rule out certain gravitational theories is the properties of gravitational waves pp|| , such as the speed of 
propagation and allowable polarization states. It was shown above how a single sphere can measure the quadrupole 
components of the strain tensor, but a scalar wave can excite both the monopolc mode and the quadrupole modes of a 
sphere 10 2]]]. By monitoring both types of modes, a single spherical detector can measure all the tensor components 
of a gravitational wave. This makes it possible for a single spherical detector to determine all of the six polarization 
states predicted by the most general symmetric metric theory of gravity p2[ . 

We can rewrite Eq. (^) in terms of the electric components of the Riemann tensor |2(| E t j = i?oioj > 



nllm 



1 

— E, 
M ' 



ij 



(16) 



where we now include both the £ = 2 quadrupole modes and the £ = monopole mode. The monopole mode of an 
elastic sphere is actually at a higher frequency than the quadrupole modes. If the source is not wide-band enough for 
detection in both of these modes, a second sphere with the monopole mode tuned to the quadrupole modes of the 
first will be needed to measure this component. If the first sphere is at relatively low frequency, one might consider 
making the second sphere hollow to keep it of a practical size p3[ |. An alternative to a second sphere is to monitor 
the n = 2 quadrupole modes and the monopole mode of a single sphere. These modes are not far in frequency from 
each other and also have relatively large cross-sections [^4|,[n| . 

Expanding Eq. ( |l6| ) into radial and angular parts we find an additional spherical amplitude ho corresponding to 
the £ = spherical harmonic. The detector response in the lab frame can now be written as 



A 



,9i 



7f 55 
92 

.94 



9o 



~9i 



.92 

73^5 

93 



9o 



9i 

93 

-^95 +.9o 



(17) 



To determine how to solve the inverse problem we need to examine the form of . It is a symmetric tensor so it 
has only six independent components. It can be written in terms of the complex Newman-Penrose parameters |2Ej | 
which allow the identification of the spin content of the metric theory responsible for the generation of the wave 



E; 



-0?* 4 - * 2 2 3*4 

3*4 5?*4 - $22 

-V83f?* 3 \/83* 3 



%/83* 3 

-6*2 



(18) 



We can divide the theories of gravity into categories using the E(2) classification scheme shown in Tab. 



_§. The 

tensor Eij is symmetric for all of these classes, thus it is orthogonally diagonalizable, but classes IIq and III5 have 
more degrees of freedom (direction plus polarization states) than we are capable of measuring with a single spherical 
detector. These two classes are often referred to as "observer-dependent" because different observers will disagree 
upon which polarization states are present. As a consequence, the polarization amplitudes for a particular observer 
must be known before the direction of the wave can be estimated. 

For the "observer- independent" classes Oo, 0\, N2, and A3, the situation is more straightforward. Oq is obviously 
uninteresting as it does not predict any gravitational waves (this class along with 0\ have essentially been ruled out by 
previous experiments |19|). We notice that E^ for the observer- independent classes can be diagonalized by a rotation 
a about the propagation axis, therefore, we can use the same arguments presented above for general relativity to solve 
for the wave direction. 

The most general observer- independent class is A3, which has 



E^ - 



-!ft*4 - $22 3*4 







3*4 5R*4-*22 















(19) 



Looking at the form of E^ we see that the same procedure for calculating the direction of the wave in general 
relativity holds for all the observer-independent classes: the eigenvector of A with eigenvalue equal to zero points at 
the source. The one exception to this statement is the case where the driving forces remain in a fixed line, for example 
3*4 = 0, 5R*4 = *22- In this situation the direction of the wave can only be determined within the plane defined by 
the two eigenvectors with eigenvalues equal to zero. 
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C. Detector response to impulsive excitations 



Impulsive excitations are often used on resonant-mass detectors to calibrate the antenna p7[ . The excitations are 
usually administered by either a short electrical burst applied to a calibrator attached to the surface or a hammer 
blow. Impulsive excitations were also used to test the analysis techniques used for experiments with the prototype 
spherical antenna at Louisiana State University (l3|Jll]. 

A radial impulse excitation can be easily described if we choose the z' axis to be along the direction of the impulse. 
By examining the quadrupole eigenfunctions of the sphere in this frame we notice that out of these modes only the 
^125 mode will be excited (other sphere modes will also be excited but their response can be removed by narrow-band 
filtering) . All of the other quadrupole modes have a vanishing radial component of their eigenfunctions at this location 
which makes their "overlap" integral with the impulse vanish. In this frame the detector response is 



A' 



M 

' 









73^ 



-i< 



(20) 



In the lab frame A is still given by Eq. (|S[> . 

Again, the direction can be found by calculating the eigenvalues and eigenvectors of the lab frame A. In the absence 
of noise, the eigenvector corresponding to the direction has a non-zero eigenvalue that is opposite in sign and twice 
as large as the two other eigenvalues. 



III. SOLUTION TO THE INVERSE PROBLEM IN THE PRESENCE OF NOISE 



We now return to the case of a gravitational wave in general relativity to solve the inverse problem in the presence of 
noise. At the end of this section we present the application of this solution to the other types of excitations mentioned 
above. For this discussion we assume that the mode channels g are independent and have normally distributed 
noise with the same variance. This is a reasonable assumption as the truncated icosahedral arrangement of identical 
transducers ideally satisfies these condition s H . In addition, several other proposals of transducer arrangements 
also produce independent mode channels |p|,^8| (but the sensitivity of each mode channel is different under normal 
conditions |2§|]). 



A. Solution for general relativity 



Noise in the mode channels g will change the eigenvalues and eigenvectors of A such that they are no longer equal 
to those of H. To gain some insight into this situation, let us consider the noise as a perturbation TV to the matrix 
H 



A = H + N . (21) 

The matrix N is constructed from the noise in each mode channel g, thus it has the same form as Eq. (j|). The 
matrix A is therefore still symmetric and traceless and has the eigenvalue equation 

The eigenvectors of H can be expanded in terms of the eigenvectors of A 

x = S ' C XX 'X , (23) 

x' 

where the matrix C is close to the identity matrix if the perturbation is small. However, since we do not know the 
values of the matrix N we cannot calculate any corrections to the matrix elements C xx i , thus the best approximation 
to x we can find is x' . 

Also from perturbation theory we see that the eigenvalue corresponding to the estimated direction of the wave is 
k x i w if the perturbation is small. Its magnitude will increase as the SNR decreases, but it should remain smaller 
than the other two eigenvalues of A for SNR > 1. Consequently the eigenvector corresponding to the estimated 



G 



direction of the wave can be selected from the three eigenvectors of A by choosing the one whose eigenvalue is 
"closest" to zero. Once x' is found, it can be used to estimate the direction of the source using Eqs. @ and (|l|). 

The perturbation approach gives us a conceptual feel for the solution, but a more rigorous proof seems necessary. 
The problem we wish to solve is to estimate the direction and polarization that makes the measured five mode channels 
g m most "look like" the expected signal from a gravitational wave, h m from Eqs. ([|). Zhou and Michelson used a 
statistical argument to justify using the least square error in their maximum likelihood method to fit for the direction 
and polarization Given the poor statistics in this estimation (only five samples) one might question their statistical 
approach, nevertheless the least squares error seems to be a reasonable choice to make under the conditions on the 
noise stated above. 



The least squares error can be written as 



= J2(9m-h m f. (24) 



The values of h+ and h x that minimize Q can be found by simultaneously solving the equations dQ/dh+ = and 
dQ/dh x = 0. Doing so using Eqs. (|J) we find 

h+ = 5i^(l + cos 2 (3) cos 27 — #277(1 + cos 2 (3) sin 27 — g 3 ^- sin 2/3 sin 7 + g^- sin 2(3 cos 7 + g$^- sin 2 f3, (25) 

£1 £i £1 £ £i 

' l x =5i cos (3 sin 27 + g 2 cos (3 cos 27 + g 3 sin (3 cos 7 + 54 sin (3 sin 7. (26) 

Note that Eqs. (^5|) and ( p6| ) are identical to Eqs. ( |l4| ) and ( |l5| ) found for the noiseless case. This connection will be 
useful below. 

We might also look for the minimum of Q with respect to the direction of the wave by taking partial derivatives 
with respect to /3 and 7. This procedure leads to very complicated non-linear equations whose solution is not easily 
obtained. For this reason we instead will look at how Q varies close to our eigenvector solution. We begin by rewriting 
Q in terms of the detector response 

3 3 

Q= 2EE(^-M 2 (27) 

j=l 1=1 



1 



2 Tr {[A- H][A-Hf) (28) 



\A-H\\ 2 . (29) 
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The inner product || A — H |j 2 can be interpreted as the distance between A and H which we know to be invariant to 
rotations. If R is the matrix that diagonalizes H such that H' = R T HR, we can write 

Q = ^\\R T AR- H'\\ 2 . (30) 

It is now clear that the least squares fit is the matrix R that minimizes the distance Q. Geometrically, this minimum 
occurs when A' = R T AR is the projection of H' onto A. Given that H' is diagonal one might guess that this 
minimum occurs when A' is also diagonal. Let us proceed to prove this conjecture. 

Let Rq be the matrix constructed from the eigenvectors of A so that RqARq = D where D is a diagonal matrix. 
Let us also assume Rq differs from R by a small rotation eW such that 

R = Ro[I + eW+^e 2 W 2 }, (31) 



where W is a skew-symmetric matrix with zeros along the diagonal. Substituting Eq. (|31| ) into Eq. (|30| ) and keeping 
terms only up to e 2 we find 

Q~h\D-H' + e[DW - WD] + l-e 2 [DW 2 - 2WDW + W 2 D}\\ 2 , (32) 
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Expanding Eq. (|32|) and remembering that the trace is invariant under cyclic permutations of the matrices in a product 
and that H'D = DH' we find 

Q ~ -\\D - H'\\ 2 + e 2 Tr{H'WDW - DH'W 2 ). (33) 

All the first order terms in e have vanished so we have proven that Q is stationary near R — R . To show this point 
is a minimum we need to evaluate the second order terms in e. 

W can be written in terms of a unit vector n representing the axis of rotation, so the square of this matrix is given 

by 

[W 2 ]ij =n i n i -%. (34) 

Recalling the procedure for deriving Eqs. ( |l4]) and ( |l5| ) in the noiseless case and that they are identical to the least 
squares minimum Eqs. ( p5| ) and ( |26| ) we can set 

H' n = -H! a = ±{D n -D 22 ). (35) 

The matrix D is also traceless so Dn + D22 + -D33 = 0. Now the second order terms can be written as 

0(e 2 ) = - (D 2 n (2 - 3n?) + D 2 2 {2 - in 2 ) + D n Z> 22 (-4 + in 2 + 3n 2 )) (36) 

Using n\ + n 2 + n\ = 1 and < n 2 < 1 we find that 0(e 2 ) is not guaranteed to be positive for all possible real 
values of Dn and D 2 2- Fortunately we may also assume we have ordered the eigenvalues of A such that D 33 is the 
eigenvalue closest to zero. Now we have an additional condition Dn = — where 0.5 < b < 2. Substituting this 
into Eq. ( |36| ) we find 

0(e 2 ) = -D% 2 (b + 1) (2b + 2 - 3bn 2 - in 2 ) . (37) 

By inspection we see that 0(e 2 ) is always positive under the conditions stated above, therefore Q is always a minimum 
near R = Rq. 

We further used a Monte Carlo type simulation to show that this point is always the global minimum of Q. For a 
wave of a given direction and polarization we calculated the spherical amplitudes h m and added a random number 
(variance a\ and zero mean) to obtain the mode channels. The direction and polarization were estimated using the 
eigenvector method as well as by numerically finding the minimum of Q from Eq. (|24|). We found the two methods 
gave identical results, even for high values of a 2 , confirming that R — Rq is a global minimum of Q. Therefore, the 
diagonal form of A is the best approximation to H' and can be used to estimate the direction and polarization of the 
wave. 

Both the maximum likelihood method and the eigenvector solution minimize the mean square error under the 
conditions on the noise stated above, therefore, produce the same answer for the estimated values. However, the 
eigenvector solution is more straightforward and computationally simple. We construct the matrix A from the mode 
channels g and compute its eigenvalues k x i and eigenvectors x' . We choose the eigenvector with eigenvalue closest to 
zero and estimate the direction of the wave using Eqs. ( |l0| ) and ([ll]). The polarization amplitudes can be estimated 
using Eqs. (p5|) and (|2^). 



B. Extensions of the eigenvector solution 

The detector response matrix A for other metric theories of gravity as well as for impulse excitations satisfy the 
symmetry arguments used in the discussion for general relativity. This means we can easily adapt the noiseless 
solutions to the case where noise is present in the same fashion. 

For observer-independent gravitational theories the method for estimating the direction of the wave is identical to 
that of general relativity: the eigenvector of the detector response matrix with eigenvalue closest to zero can be used 
to estimate the direction of the source. Observer-dependent theories require prior knowledge of the polarization states 
of the wave before any estimate of the direction can be made. Once these are known it should be straightforward 
to adapt the eigenvector technique to estimate the direction. In the case of an impulsive excitation, the eigenvector 
corresponding to the direction has an eigenvalue that is opposite in sign and greater in magnitude than the two other 
eigenvalues. Converting the eigenvectors to a direction again comes from Eqs. © and |H]). 
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IV. DISCUSSION 



The eigenvector solution is very convenient in that the inverse problem is reduced to solving a trivial eigenvalue 
problem. The solution is computationally simple, making this technique very efficient for use in an automated data 
analysis system. This feature may be important if one considers using a large number of candidate gravitational wave 
events in a coincidence exchange between several detectors where the source direction is used as a criterion to veto 
excess coincidences. 

The main restrictions on the eigenvector solution is that the mode channels must be independent and the noise 
normally distributed with equal variance. These restrictions can ideally be satisfied for a number of transducer 
arrangements We found that the eigenvector solution corresponds exactly to the maximum likelihood method 

under these conditions. It may be possible to apply this solution when the noise is not gaussian or is different for 
each mode channel, however further research is necessary to verify this extension. 

Through a number of numerical simulations as well as examination of the work of others [||J2^] we found that the 
errors due to the noise on a direction estimation are independent of the source location and wave amplitude for a given 
SNR. However, the estimation of the polarization amplitudes using Eqs. (^5|) and (^) lead to direction dependent 
uncertainties. For example, Fig. |^ shows the variance on the polarization angle a = tan -1 (/i x //i + ) for a range of SNR 
and several values of f3 found from a Monte Carlo type simulation. Notice that the variance increases for low values 
of [3. This realization is disturbing given that a spherical antenna is equally sensitive to waves from all directions and 
polarizations. 

One might consider using a different coordinate system to try to avoid the directions with very poor estimates of 
the two polarization amplitudes. For example, use the xyz-convention of the Euler angles where the first and the 
last rotations are not the same. This actually will not solve our problem, but instead change the directions in the 
sky which lead to the poor estimates. If we transform back from this coordinate system to the y-convention we just 
reintroduce the errors and thus have gained nothing. 

This dependency on the source direction is not unique to a sphere, a network of bars or interferometers will also 
suffer from this problem |l2| . This leads us to believe that we are excluding a piece of information from our procedures. 
The solution may lie in using the information from the two other eigenvectors of the detector response. In the above 
derivations these eigenvectors were simply discarded, but they also contain information about the gravitational wave 
that may eliminate these direction dependent errors. This approach will be the topic of a future paper p9fl . 

While there are a few limitations to the eigenvector solution, its simplicity makes it easily extendable to other 
types of excitations. As discussed above, impulsive excitations can be located using this technique. As a practical 
example we recall that this solution was successfully tested on experiments with the LSU prototype spherical antenna 
[]14j . This practical confirmation of its validity gives us the confidence that it can be implemented on a real spherical 
antenna searching for gravitational waves. 
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FIG. 1. The y-convention of the Euler angles. 
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FIG. 2. The results of a numerical simulation describing the variance of the polarization angle a for a range of SNR and 
several values of the direction angle f3. Each line was computed by a 500 trial Monte Carlo simulation for 100 logarithmically 
spaced SNR for the corresponding value of (3. 



Class Allowable polarization states Example 

1% $22, *4, *3, *2 Most general 

Ills $22, *4, *3 Kaluza-Klein 

A^3 $22, *4 Brans-Dicke 

N2 *4 General relativity 

Oi $22 Purely scalar 

Oo None No wave 



TABLE I. The E(2) classification scheme. 
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